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Abstract. Blind 21cm surveys in the local universe have shown that 
the local H I mass density, J)hi> is dominated by luminous, high surface 
brightness, spiral galaxies. On the other hand, surveys for host galaxies of 
damped Ly-a systems have not always been successful in finding bright 
spiral galaxies. From an analysis of 21cm aperture synthesis maps of 
nearby galaxies we show that this apparent contradiction can be resolved 
by realizing that the H I mass density is dominated by L* galaxies, but the 
H I cross section near the DLa threshold is more evenly distributed over 
galaxies with a large range in luminosity, gas mass, and surface bright- 
ness. The distributions of column densities and impact parameters of 
optically identified and non-identified DLa host galaxies in the literature 
and the H I maps are qualitatively in agreement. Due to poor number 
statistics of low redshift DLa systems, there is no firm indication that the 
redshift number count of low redshift DLa systems is inconsistent with 
that calculated from the nearby galaxy population. 



1. Introduction 

Blind 21cm surveys and studies of damped Ly-a (DLa) systems use completely 
unrelated observational techniques but they both can be employed to address 
similar scientific questions: what is the cosmological mass density of H I (Ohi), 
what sort of galaxies contribute most to Ohi> and how does Ohi at z = relate 
to that at high zl The main advantage of the 21cm surveys is that the total H I 
mass of the detected cloud or galaxy can be readily derived, and after follow- 
up with synthesis instruments a 3-dimensional data cube can be obtained from 
which a detailed map of the H I distribution and a velocity field can be extracted. 
Unfortunately, the neutral hydrogen 21cm emission line is weak, and blind 21cm 
surveys are limited to the very local {z < 0.1) universe. In contrast, distance is 
not a limiting factor in the identification of DLa systems, since their detection 
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depends on the brightness of the background source against which the absorber 
is observed. However, the nature of these systems is more difficult to determine 
because the absorption spectrum gives information along one very thin sight-line 
through the absorber, and the total H I mass of individual absorption systems 
can not be determined. Furthermore, because DLa systems are often found 
at high redshift, determining their optical properties requires extensive, time- 
consuming imaging programs. 

One of the major challenges of 21cm and QSO absorption line research is 
to reconcile the results from blind H I surveys and DLa surveys at low redshift. 
The aim of this paper is to briefly review the results from both fields, and test 
whether the results can be brought in agreement. 



2. Results from 21cm and DLa surveys 

The most extensive blind 21cm surveys to date are the Arecibo H I Strip Survey 
(AHISS, Zwaan et al. 1997), the Arecibo Dual Beam Survey (ADBS, Rosenberg 
& Schneider 2000), and the HI Parkes Sky Survey (HIPASS, Staveley-Smith et al. 
1996, Kilborn 2001). The main results from the AHISS, relevant to the present 
discussion, are that (i) the local H I mass density is dominated by H I rich, high 
surface brightness, luminous spiral galaxies (Zwaan et al. 1997, Zwaan, Briggs 
& Sprayberry 2001). Galaxies with H I masses below 1O 8 ' 5 M make up only 
about 20% of the neutral gas density, and low surface brightness (LSB) galaxies 
with central 5-band surface brightness > 23.0magarcsec -2 contribute a similar 
fraction, and (ii) where there is H I there are stars. No convincing instance of 
star-less extragalactic H I cloud, not gravitationally bound to a parent galaxy, 
has yet been identified. The results from the ADBS and preliminary results from 
HIPASS are in agreement with these conclusions. 

Thus, if most of the H I in the local universe is in these conspicuous galaxies, 
and at z > DLa systems contain most of the H I, this would naively lead to 
the conclusion that DLa host galaxies should be predominantly bright spiral 
galaxies and thus easily identifiable. However, this does not appear to be the 
case. The tally so far is that optical imaging of z < 1 DLa systems has turned 
up 5 LSB irregulars, 3 compact dwarfs, and ±5 'normal' spiral galaxies (Le Brun 
et al. 2000, Nestor et al. 2001). In ±2 cases (depending on which reference to 
choose) a host galaxy has not been identified. Ha emission has been searched 
for in only one DLa system by Bouche et al. (2000). They did not detect Ha 
emission in their deep HST imaging of a z = 0.656 DLa system, which allows 
them to put an upper limit of 1.9 Mq/yt on the global star formation rate in 
DLa system, which is lower than the typical value measured in L* field spiral 
galaxies. A third method of studying properties of DLa systems is to search for 
21cm emission in order to measure the total cold gas content of the DLa system. 
The problem with these surveys is that the 21cm H I emission line is extremely 
weak so that with present technology these surveys are limited to the very local 
(z < 0.2) Universe. Two attempts to detect H I in emission in DLa systems have 
been made by Kanekar et al. (2001) and by Lane et al. (2000). Both searches 
resulted in non-detections, which allowed the authors to put upper limits to the 
H I masses of approximately Mm = 3 x 10 9 M Q (Hq = 65kms~ 1 Mpc -1 ), which 
is one third of the H I mass of an L* galaxy. A reverse strategy was followed by 
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Bowen et al. (2001), who identified a DLa absorber in a nearby (z = 0.01) very 
gas-rich LSB galaxy. 

The conclusion from these surveys is that low redshift DLa host galaxies 
display a wide range in galaxy parameters. This is a conclusion very similar to 
that from Nestor et al. (2001) who discuss the results of an optical and NIR 
imaging survey of five QSO fields with identified DLa systems with z < 0.53. 
Obviously, it are not just the bright, high surface brightness spiral galaxies that 
dominate the cross section for TVhi > 2 x 10 20 cm -2 . Is this a contradiction to 
the low redshift 21cm emission line survey results? To answer this question it is 
essential to study the Nm distribution in nearby galaxies in more detail. 

3. Measuring the local H I column density distribution function 

Determination of the local H I column density distribution function requires 
21cm maps of a fair sample of galaxies that is representative of the nearby galaxy 
population. Few such samples exist, because traditionally radio astronomers 
have chosen to study galaxies with extended H I disks, so that the rotation curve 
can be measured out to large radii, or galaxies with peculiarities or participating 
in interactions. The currently available data from H I selected galaxy samples 
are not appropriate because the spatial resolution of the survey instruments 
(Arecibo or Parkes) is too coarse to make detailed H I maps. 21cm follow-up 
on some of these samples (AHISS, ADBS) has been done with the VLA in the 
most compact (D) configuration, which only barely resolves the galaxies. Using 
these low resolution H I maps at a particular limiting column density contour 
leads to an overestimation of the effective cross-section for DLa absorption. 
The exception is the ongoing work of Ryan- Weber et al. (2001), who use the 
Australian Telescope Compact Array (ATCA) to study a large sample of nearby 
galaxies taken from the HIPASS catalogue. Eventually, this sample will yield a 
valuable measurement of the local H I column density distribution function. 

In Zwaan, Verheijen & Briggs (1999) we presented a measurement of the 
z = H I column density distribution function [f(N)], derived from deep WSRT 
21cm maps of a volume-limited sample of galaxies in the Ursa Major (UMa) 
cluster (Verheijen & Sancisi 2001). We argued that this sample of galaxies is 
representative of the field because (i) the collapse time of the Ursa Major cluster 
is close to a Hubble time which implies that few interactions have taken place, 
(it) the H I maps show no signs of H I depletion, and (in) the morphological 
mix of the galaxy population in the Ursa Major cluster is similar to that in the 
field. The sample is very suitable for calculating f(N) because all galaxies are 
at approximately the same distance, which yield a uniform spatial resolution of 
~ 1 kpc. The sample consists of 49 galaxies, and is complete to an absolute 
magnitude Mb = —16.9, roughly twice the luminosity of the Small Magellanic 
Cloud. 

The UMa f(N), presented in Figure 2 of Zwaan et al. (1999), was a consid- 
erable improvement over the function previously derived by Rao &: Briggs (1993), 
who used Gaussian fits to the radial H I distribution of 27 galaxies with large 
optical diameters, as measured with the Arecibo Telescope. The UMa analysis 
was based on much higher spatial resolution observations (~ 15" vs. 3'), the 
sample contained more galaxies, and no modelling was applied. The UMa f(N) 
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log N HI > 19.8 log N H , > 20.3 log N HI > 20.8 log N HI > 21.3 

Figure 1. Cross-sectional area of HI in the UMa galaxy sample as 
a function of 5-band absolute magnitude. The different colors corre- 
spond to different column density cut-offs, where the fat, red line is for 
column densities logiVHi > 20.3, which is the limit above which an HI 
absorbing system is called a damped Ly-a system. The three panels 
show the effect of using different assumed shapes of the galaxy lumi- 
nosity function. The numbers show the values of 100 x dN/dz, where 
dN/dz is integrated over the the whole range in Mb over which the 
function is evaluated. Errorbars are calculated using Poisson statistics. 



follows a A^hi 1 distribution for Nyj < 10 21 cm 2 , and a Nm 3 distribution for 
N m > 10 21 cm" 2 . 



4. The distribution of cross-sectional area 

Here we re-analyze the Ursa Major sample to make a more specific compari- 
son between the available DLa data and the information from nearby galaxy 
population. In Figure 1 we show the 'volume density of cross-sectional area' 
for different column density cut-offs. This cross-sectional area per volume is 
calculated by multiplying A(Mb) x <&(Mb), where A{Mb) is the area of the 
H I above a certain column density limit, as a function of absolute magnitude 
Mb, and $(Mg) is the galaxy luminosity function. The result of this product 
is dN/dz x Hq/c, where dN/dz, the number of absorbers per unit redshift, is a 
familiar quantity in QSO absorption line studies. To account for the fact that 
the sample excludes galaxies with inclinations i < 45°, we have de-projected 
the H I maps, and recalculated the area function A(Mb) for the sample as if 
it were randomly oriented with inclination i < 45°. These values have been 
incorporated in the total value of A(Mb)- 

To illustrate the effects of different shapes of the luminosity function, we 
show the distribution of A(Mb) X $(Mg) for three different values of the faint- 
end slope of the Schechter function: a = —1.0, a = —1.2 and a = —1.5, as 
indicated by the miniature luminosity functions in the upper left corner of each 
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Figure 2. Cross-sectional area of HI in the UMa galaxy sample as 
a function of H I mass. (Curves are labeled as described in Fig. 1.) 
The H I mass function used for the normalisations are illustrated in 
the upper left corner of each panel. 



panel. The value of a = —1.2 is closest to the recently published results from 
the 2dF survey (Folkes et al. 1999) and the Sloan Digital Sky Survey (Blanton 
et al. 2001). 

The main point here is that the distribution of cross section for log A^hi > 
20.3 is much less peaked around L* galaxies than the distribution of H I mass 
density (Figure 3 from Zwaan, Briggs, &: Sprayberry 2001). Indeed, the H I 
cross section in galaxies with Mb = —16 is only slightly lower than that in 
galaxies with Mb = — 20. The consequence of this observation is that if an H I 
column density > 10 20 ' 3 cm -2 is encountered somewhere in the local universe, 
the probability that this gas is associated with a L*/40 galaxy is almost as 
high as for association with a L* galaxy. On the other hand, the highest column 
densities (> 10 2L3 cm -2 ) are largely associated with the most luminous galaxies. 
Small cross sections of very high column density gas can still account for high H I 
masses, which explains why the total H I density is dominated by L* galaxies. 

A similar plot can be constructed for the cross-sectional function 
of H I mass. This is shown in Figure 2, where we plot A(Mb) x 0(Mhi), 
where 0(Mhi) is the H I mass function. Again, the results for three different 
values of the faint-end slope are shown. The H I mass fuction based on the 
Arecibo H I Strip Survey (AHISS) discussed by Zwaan et al. (1997), has a 
slope a = —1.2, corresponding to the middle panel. Steeper H I mass functions 
have been advocated by Rosenberg & Schneider (2001) on the basis of the ADB 
Survey. 

The dependence of interception cross section on HI mass of the host in Fig. 2 
shows a similar behavior to that of optical luminosity (as in Fig. 1), although 
less pronounced. Most of the cross section of DLa column densities is associated 
with large H I mass galaxies, but the contribution from lower mass galaxies is not 
negligible: galaxies with Afni < 10 9 M Q contribute ~ 30% to the cross section. 
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Figure 3. Same as Figure 1, but now as a function of i^-band central 
surface brightness. The same luminosity function as in Figure 1 have 
been used to calculate normalisations. The arrow with 'Freeman' indi- 
cates the -RT-band equivalent of the 'Freeman- value', the mean central 
surface brightness of high surface brightness spiral galaxies. 



Finally, Figure 3 shows the cross-sectional area as a function of K-band 
central surface brightness. Weights to the individual galaxies have been calcu- 
lated using optical luminosity functions, similar to those in Figure 1. This is a 
more robust method of determining the normalisation than using the central sur- 
face brightness distribution function (the galaxy density as a function of central 
surface brightness), which has not been measured to great accuracy. Concentrat- 
ing on the central panel, it appears that galaxies two magnitudes dimmer that 
the canonical 'Freeman- value' (Freeman 1970) contain a comparable amount of 
cross section at DLa column densities to those with the Freeman central surface 
brightness. Again, the distribution function for the highest column densities is 
more clearly peaked around high surface brightness galaxies, which explains the 
dominance of these galaxies in Ohi- 

The conclusion from these three figures is that the probability distribution 
of cross-sectional area of DLa column densities as a function of various galaxy 
parameters is much less peaked at L*-type properties than the H I density. This 
explains why surveys for DLa host galaxies turn up a population of galaxies 
with a large spread in luminosity, morphological type, and surface brightness. 
Rosenberg & Schneider, and Turnshek et al. have presented similar plots of 
cross section vs. luminosity at this meeting. 



5. Redshift number density 

The above analysis shows that 21cm surveys and DLa results are in qualitative 
agreement. The remaining question is whether dN/dz, the number of absorbers 
per unit redshift, is quantitatively consistent between the two techniques. 

Rao & Turnshek (2000) recently undertook a HST survey of low redshift 
(z < 1.65) Mgn selected systems, with the aim of measuring the incidence of 
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DLa systems. The underlying assumption that makes this statistical approach 
work, relies on the empirical fact that all DLa systems are Mgn absorbers. 
Rao & Turnshek conclude that dN/ dz shows no indication of evolution over the 
redshift range 0.5 < z < 4, but the z = point measured from 21cm surveys 
is significantly lower. However, Briggs (2001) argued that the Mgn statistics 
used to bootstrap the DLa dN/dz values, should be limited to those systems 
with W/" A2796 > LOA, since the Mgn-selected DLa systems are drawn almost 
entirely from the large equivalent width Mgn sub-sample. Because Steidel & 
Sargent (1992) found that at low redshift these large Wq 2796 systems evolve 
significantly faster than the low Wq 2796 systems, the lowest z Rao & Turnshek 
values of dN/dz would be lower, which could lessen the discrepancy with the 
z = point substantially. Interestingly, Churchill (2001) performed an unbiased 
survey for low-z Mgn systems and used the same statistics for Mgll systems as 
Rao & Turnshek did, and found that dN/dz for DLa absorption is 0.08^q'q5 
at (z) = 0.05. This value agrees very well with a no-evolution extrapolation of 
the Rao & Turnshek results. Unfortunately, small number statistics (Churchill's 
result is based on four systems) causes large uncertainties on the derived values, 
which makes the interpretation very difficult. The f2ni studies that depend 
on intermediary Mgll lines contrast sharply with lower values from the more 
traditional DLa survey methods (Jannuzi et al. 1998, Lanzetta et al. 1995), 
although they also suffer from small number statistics. 

Using the UMa sample and the H I mass function from Zwaan et al. (1997) 
with a faint-end slope of a = —1.2, we find for DLa column densities a value of 
dN/dz(z = 0) = 0.042 ± 0.015. The quoted uncertainty on this value is based 
on Poisson statistics, and does not take into account the uncertainties on the 
H I mass function parameters. If we use the isophotal S-band galaxy luminosity 
function parameters from Blanton et al. (2001), and account for the galaxies in 
the UMa sample not detected in H I, we arrive at the same value for dN/dz. This 
value is higher than the one found by Rao & Briggs (1993) (dN/dz = 0.015), 
probably because these authors limited their calculation to large, optically bright 
galaxies, whereas our calculation is based on a volume-limited galaxy sample, 
which also includes dwarf and LSB systems. 

At present it is unclear whether the DLa redshift number densities at z ~ 
derived from 21cm emission maps and Mgn selected systems agree. The error- 
bars on the Churchill value are sufficiently large to make them in agreement 
with our estimate (less than la difference). Obviously, much better statistics 
are required to appraise the results from the different methods. Rosenberg & 
Schneider (2001) use their VLA D-array imaging of H I selected galaxies to find 
dN/dz = 0.07 ± 0.01. There are two reasons why this value is higher than the 
one we find. Firstly, the H I mass function that these authors use has a faint-end 
slope of a = —1.5, whereas we use a = —1.2. This steeper H I mass function 
causes low mass galaxies to contribute more cross section. Secondly, their anal- 
ysis is based on lower resolution data of galaxies at larger distances. This might 
lead to an overestimation of the cross sections. 
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Figure 4. The two-dimensional probability distribution of cross sec- 
tion in the (Nyu,R) plane, where R is distance from the centre of the 
galaxy. The contours are from the H I maps from the UMa sample, and 
the dots represent data from DLa surveys. The shaded area indicates 
the region that corresponds to DLa column densities. 

6. Comparing properties of DLa systems and nearby galaxies 

In Figure 4 we show the 2-dimensional probability density distribution of cross 
section in the (iVm, R) plane, where R is distance from the centre of the galaxy. 
We have used the UMa sample to calculate the cross-sectional area contributed 
by each element dNmdR, on a fine grid in the (Nm,R) plane. To increase 
the signal-to-noise in the figure, we smoothed this probability distribution with 
a Gaussian filter with a = 0.1 in the logiVHi direction and a = 1 in the R 
direction. An H I mass function with a faint-end slope a = —1.2 was again used 
to calculate the weights of the individual galaxies. The contour levels are chosen 
at 10, 30, 50, 70, and 90% of the maximum value. 

There are a number of interesting features in this plot. First, it shows that 
the largest concentration of H I cross sections in galaxies in the local universe is in 
column densities 20.6 < log A^hi < 21.1 and galactocentric radii 3 < R < lOkpc. 
The observational fact that the H I distribution in galaxies often shows a central 
depression can be seen in this figure by the compression of contours at R ~ 0. 
Not surprisingly, the outer contour shows that high column densities are mostly 
seen at small R, but it also shows that column densities lower than log Nm = 20.6 
are rarely seen in the extreme centres of galaxies. (This might depend on the 
resolution) . 

In order to make a direct comparison between 21cm results and DLa pa- 
rameters, we plotted on top of the contours the data points from DLa surveys 
taken from Le Brun et al. (2000) and Nestor et al. (2001). Note that these are 
the lowest redshift DLa systems known, but the median redshift of the com- 
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bined sample is (z) = 0.5, significantly higher than the median z = 0.003 for the 
UMa sample. Evolutionary effects might therefore be important, but at present 
a lower redshift DLa sample is not available. 

Although the statistics is still poor, there seems to be an agreement be- 
tween the distribution of the DLa points and the contours from the 21cm maps. 
When more measurements of impact parameters and column densities of low- z 
DLa systems become available, this method will provide a powerful tool to test 
whether the properties of DLa systems and nearby galaxies are in agreement. 



7. Concluding remarks 

We have shown that the properties of DLa galaxies are in good agreement with 
the expectations based on 21cm images of nearby galaxies. Although most of the 
H I mass density is locked up in luminous spiral galaxies, the H I cross section is 
much more spread out over different galaxy properties. Also the distribution of 
column densities and impact parameters are in general agreement between DLa 
systems and local galaxies, although the statistics is still poor. 

There are a number of caveats that should be kept in mind when making 
comparisons between QSO absorption line systems and galaxies. First, the spa- 
tial resolution of the aperture synthesis instruments used to derive H I column 
density maps of nearby galaxies might have a influence of the derived column 
density distribution functions. Beautiful examples of fine-structure in the H I 
distribution in galaxies can be found, for example, in the maps from the Southern 
Galactic Plane Survey (McClure-Griffiths et al. 2001), in the Small Magellanic 
Cloud maps of Stanimirovic et al. (1999), or in Braun's (1995) high resolution 
imaging of nearby spiral galaxies. The results presented here are based on much 
coarser resolution, which means that most of this fine structure is averaged out. 

A second point of concern is the small number of appropriately studied 
galaxies in the available samples. The cross section distribution functions are 
derived from 49 galaxies, which implies that there are on average 10 galaxies per 
bin in Mb, Mm, or surface brightness. The calculation of cross-sectional areas 
implicitly assumes that within each bin the galaxies are randomly oriented, but 
for these small samples that might not always be justified. Also, it would be in- 
teresting to calculate the functions for a larger range of galaxy parameters. Both 
these issues will be addressed in a future paper in which we present an analy- 
sis of the column density distribution function of a much larger galaxy sample 
taken from the WHISP survey (Westerbork observations of neutral Hydrogen in 
Irregular and SPiral galaxies). 

Finally, it is important to keep in mind the biases that are introduced in 
assembling DLa galaxy samples. As Turnshek et al. (2001) stress, the identifica- 
tion of a galaxy with the DLa absorber is more probable if the impact parameter 
between the sight-line to the background QSO and the galaxy is small. On the 
other hand, if the impact parameter is very small, the identification of a galaxy 
might be hindered by the dominating QSO light. DLa samples might be biased 
due to the effect of dust obscuration in the DLa host galaxies. Pei et al. (1999) 
argue that these 'dusty sight-lines' might drop out of the QSO sample that is 
used to identify the absorbing systems. 
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